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Abstract. For the intermediate-type barred galaxy NGC 2336 
stationary models are constructed which reproduce in a consis- 
tent manner the observed distribution of the luminous matter 
and the observed gas kinematics in those regions affected by 
the bar We present 2D fits to the observed iV/i?-band lumi- 
nosity distribution that consist of three components; a bulge, a 
bar, and a disk. The brightness distribution of each component 
is converted into an underlying mass distribution by means of 
a suitable A//i-conversion. The resulting coadded potential of 
NGC 2336 is implemented into a numerical code for the com- 
putation of closed orbits for gas clouds (Hll-gas). Using the 
resulting tracks, the phase space accessible to the models is 
examined with regard to the main orbit families. For different 
orbit energies complete sets of closed orbits are computed. By 
projection to the reference frame of the galaxy, artificial rota- 
tion curves for every model are obtained and are compared with 
the observed rotation curves of the Hll-gas. 

In an iterative procedure, the parameters of the NGC 2336- 
models are optimized by computing and evaluating a large 
number of parameters. The result is a final model that repro- 
duces the morphological structure of NGC 2336 as well as the 
observed kinematics of the Hll-gas. The parameter values from 
the morphological decomposition and those needed to fit the 
Hll-rotation curves best are in exellent agreement. The effects 
of changing single parameter values and possible error sources 
are discussed in detail. It turns out that the kinematics of the 
warm Hll-gas of NGC 2336 can be explained without consid- 
ering hydrodynamic effects, even in the central regions. 

Key words: galaxies: kinematics and dynamics - galaxies: 
photometry - galaxies: structure - stars: kinematics - ISM: 
kinematics and dynamics 



1. Introduction 

Due to their strong non-axisymmetric potential, barred galax- 
ies exhibit peculiar stellar and especially pecuUar gas kine- 
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matics when compared to normal disk galaxies without bar 
Since deconvolving their non-axisymmetric velocity field is 
difficult, the computation of consistent multi-component sta- 
tionary models for these objects has proven to be a suitable 
method to gain further insight into the internal mass distribu- 
tion and kinematics. This stationary modelling procedure can 
answer the question whether the observed kinematics of the 
galaxy agree with the model predictions derived from the ob- 
served distribution of luminous matter (assumed that plausible 
corrections for the presence of dark matter can be appUed). 

For the following studies the northern hemisphere object 
NGC 2336 was chosen. NGC 2336 is an intermediate-type spi- 
ral of Hubble-type SAB{r)bc with a prominent bar. This 
galaxy has not been subject to extensive morphologic al and 
kinematical studies yet. Basic parameters given by Tully ( 1988| ) 
are listed in Tablejlj. NGC 2336 contains sufficient amounts of 
warm Hll-gas for emission-line-spectroscopy. As is also re- 
quired for kinematic studies and deprojection, the bar lies in 
a non-degenerate position with respect to line of sight, i.e. its 
apparent major axis does not coincide with the kinematic line 
of nodes (denoted LON hereafter). The inclination of i = 59° 
ensures a fairly precise determination of the apparent axes ra- 
tios of disk, bulge and bar. 

Former photometric studies of NGC 2336 w ere ma inly re- 
stricted to the HI observations of van Moorsel (1983). In HI, 



this galaxy reveals a very regular morphological structure de- 
void of any kind of anomalies or major asymmetries, except for 
the nearly complete lack of HI in the central regions. The HI 
distribution supports the image of a mainly undisturbed spiral 
structure in the outer part of the disk with numerous star form- 
ing regions. Although NGC 2336 belongs to an apparent pair of 
galaxies (together with IC 467) with a projected linear distance 
of 135kpc, its undisturbed disk does not exhibit any distinct 
sign of recent interactions. No morphological structures con- 
necting these objects have been found. Due to the large beam 
size used for the Hl-observations {d — 35"), the inner part of 
the disk with the bar is not resolved. The HI disk isophotes at 
large radii remain undisturbed down to the detection limit. With 
lower threshold values of uhi > 1-2 • 10'^"atoms/cm-^, the disk 
extends up tow 55kpc (8.2') from the center Recent studies 



by Martin ( |1995| ) based on POSS plates yielded the basic mor- 
phological parameters for the bar listed in Table y. 
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Table 1. Basic parameters of NGC 2336 as given in Tully 
(1988). The quantities are defined as follows: a2000j /32000-' 
right ascension and declination at the epoch 2000.0, i: in- 



jb,i 



clination from face-on, Br^ : blue apparent magnitude cor- 
rected for reddening due to internal and external absorption. 



I?25 ■' diameter at the 2^mag/arcsec^ blue isophote, adjusted 
for projection and obscuration effects, R: distance based on 
Hq = lOOkm/sec, Lb: intrinsic blue luminosity, Mt: total 
mass of NGC 2336, Mhi/Mt: ratio of HL -mass to blue lu- 
minosity, Mt/Lb ratio of total mass to blue luminosity. 



According to the results of the Hl-studies by van Moorsel 



(1983), NGC 2336 does not exhibit any kinematic peculiari- 
ties up to a distance of 5' (33kpc) from the center, i.e., the HI- 
velocity field is typical for a disk-dominated galaxy. A distor- 
tion of the velocity field at 40" north to the center is caused by 
a sudden drop of the Hl-column density and does not indicate 
non-circular motions in that region. The velocity field shows 
that the southern part of the galaxy is inclined towards us. 
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0.70 
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0.09 
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0.59 
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0.17 



Table 2. Morphological parameters for NGC 2336. The re - 
sults are based on P OSS-plate measurements by Martin ( 1995 ). 
Parameters are: PAdisk: position angle of the apparent disk 
major axis, atar, bbar: major and minor axis of the bar, 
b/dbar, Lf,: apparent bar axis ratio and relative bar length 
compared to I?25' b/a{i), Lb{i): the same as above, but ad- 
justed for deprojection. 



Because of the low resolution and the lack of HI just in the 
inner region, the HI velocity field is not suited to study the per- 
turbation of the velocity field due to the bar potential. Optical 
spectrograms which show the kinematics of stars and (warm) 
gas with a much better resolution have not been available up to 
now. 

The aim of this paper is the quantitative understanding and 
modelling of the bar-perturbed velocity field in NGC 2336. For 
that purpose we study the NIR morphology of this galaxy, rep- 
resenting the ditribution of the luminous matter The total po- 
tential of disk, bar, and bulge is calculated. We consider the 



Hll-gas to be in a stationary motion in this potential. The sub- 
sequent kinematical modelling leads to artificial rotation curves 
which are compared with the observed velocity field. The op- 
timal parameters for disk, bar, and bulge, as well as the intrin- 
sic geometry of these components, are obtained by an iterative 
study of many different models and their closed orbits. 

The paper is divided into the following parts: Section 2 ad- 
dresses the photometric observations, section 3 introduces the 
procedure for constructing multi-component 2D mass models 
and their subsequent deprojection. Section 4 presents the data 
reduction process and the detailed results of the spectroscopic 
observations of both, the gaseous and the stellar component 
of NGC 2336. A basic outline of the numerical representation 
of the potentials used for the orbit integrator (which computes 
closed orbits by means of a numerical FORTRAN-/C-program) 
is given in Section 5. Section 6 refers to the different Lind- 
blad resonances in the NGC 2336-models, to the different orbit 
families and their relative contribution to the total number of 
orbits in a model. Model units are discussed in Section 7, while 
the basic instrument for phase space analysis - the Poincare- 
Surfaces of Section - are introduced in Section 8 where the 
characteristic features of the NGC 2336-phase space are dis- 
cussed. Finally, section 9 deals with the variation of indepen- 
dent model parameters which lead to the optimal model solu- 
tion in Section 10. Results are discussed in Section 11. 

2. Surface Photometry 

2.1. NIR Images 

To avoid wrong estimates for the relative contribution of the 
model components of NGC 2336 to the total luminosity, all 
models are based on NIR images (J-band) where dust ab- 
sorption is much less important than in the optical wavelength 
range. Additionally, galactic structures appear much more reg- 
ular due to the smoother distribution of cool giant stars whose 
light is mainly traced in the NIR (for a review see Frogel et 
al. 



1996) 



The images used for our studies were obtained using the 
M AG/C-camera at the MPIA 2.2m-telescope at Calar Alto 
(Spain) during September 1995. MAGIC consists of a 256 x 
256 NICMOS-chip and provides a field of view of 171" with a 
resolution of 0.67"/pixel when mounted at the 2.2r7T,-telescope 
(//8). Images were obtained in the J- (1.2/im) and A'-band 
(2.2/im). 

By repeating a routine which centers the telescope alter- 
nately on the object and neighbouring sky fields, 48 object ex- 
posures and skyflats with exposure times of t^xp = lOsec were 
obtained, resulting in a total exposure time of t^xp = 480sec. 
For the skyflats, the telescope was moved 25' away from the 
object field in several directions. An additional telescope offset 
of a few pixels between the repetitions of the routine was also 
applied to correct for bad pixels on the detector chip. Several 
domeflat series were exposed to correct for illumination effects 
on the NICMOS array. 
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Fig. 1. J- f/e/fj one/ K-band image (right) of NGC 2336 with a resolution of Q.Ql" /pixel. The disk extends well beyond the area 
of the NICMOS-detector Due to the stronger airglow emission, the S/N-ratio in K is lower than in J. 



2.2. Data Reduction 

The single exposures of NGC 2336 were reduced by a standard 
reduction process: A 2D sky pattern is constructed by using 
a vertical median filter for the sky fields that were obtained in 
different offset directions. The resulting image is normalized to 
zero. The sky contribution is subtracted by linear interpolation 
of the sky level between the sky exposure before and after each 
object frame. By median-filtering the domeflat series, a master 
domeflat is created and normalized to one. The 2D sky pattern 
is subtracted from the object frame and the result is divided by 
the master domeflat. The pre-reduced object frames are then 
centered on field stars and then co-added by a K-cr-clipping al- 
gorithm to avoid the disadvantages of a simple addition (bad 
pixels remain) or a simple median filtering (quantitative distor- 
tion of the images, lower 5/A^-ratio). The resulting images of 
NGC 2336 are shown in Fig. ^ 



zero models for disk and bar, only the bulge model will make 
use of a 3D profile. 

The equation for the flux model of the 3 components is 



Ftot{x,y) = Fbulge{x,y) + Fdisk{x,y) + Fbar{x,y). 



(1) 



3. Modelling the Mass Distribution 

The aim of the mass modelling is the determination of the un- 
derlying density distribution from the observed light profile of 
NGC 2336. The procedure works in several steps: First we de- 
termine the light distribution from the NIR images and decom- 
pose it into the three components disk, bar, and bulge. By it- 
eration the optimal parameters of each component like scale 
lengths, central flux density etc., were determined by fitting an 
appropriate model. Since we used NIR images, this morpho- 
logical model gives a good representation of the st ellar light 
and also mass distribution (see Frogel et al., 1996| ). In order 
to keep the number of free parameters small we use thickness Alt, 



For the derivation the mass distribution, one has to take into ac- 
count individual mass-to-light ratios M/L for the components: 

Mtot — Mbulge + Mdisk + ^har (2) 

= Fbulge * {M/L)bulge + Fdtsk * {M / L)disk (3) 

+ Fbar * {M/L)bar 

For this step we have to consider the dark matter halo compo- 
nent which has an increasing dynamical influence with increas- 
ing radius. However, for spiral galaxies of medium and high 
luminosity the dynamical mass within the optical diameter of 
the galaxies is strongly dominated by the disk (e.g. Salucci et 
al, 1991). In the inner zone the dark matter does not play a 
dominant role. Since we are mainly interested in the dynamical 
influence of the bar, we refrain from modelling the dark mat- 
ter component explicitly and include its contribution in (M/L) 
of the disk (the bar). This means that our final {M/L)disk 
{{M / L)bar) should not be directly compared with observed 
[M/L) values for spiral disks (bars) since, besides the pop- 
ulation effect, it contains a contribution from the dark matter 
halo. This procedure was chosen to include the dark matter in 
the corresponding potentials in a simple way, without making 
further assumptions about the spatial distribution of the dark 
matter Therefore, we hereafter denote the mass-to-Ught ratios 
by Cbuige,Cdisk, Cbar- Equation 3 becomes 



F 



*Cb 



bulge * ^bulge 



Fdisk * Cdi. 



(4) 
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Fig. 2. Left: J -band image of NGC 2336 (0.67"/pixel), center: 2D three-component model of the observed luminosity distribution, 
right: residuals remaining after subtraction of the model from the original image. 



Cbuige is nearly unaffected by dark matter contributions. Since 
the morphological model is based on J-band images, we set 



{M/L)j^l,ulge — Cbuige — 1- 



(5) 



In a next step we define Cdisk and Cbai- as relative mass-to-light 
ratios, compared to the bulge. We start the dynamical models 
with a set of start values for Cbuige, Cdisk, Char- These pa- 
rameters will be optimized by an iterative comparison between 
the kinematical models and the spectroscopic observations, see 
section 9. The calibration of the total mass will be performed 
via the observed circular velocities in the outer parts of the disk, 
see section 10. 



3.1. Profile types 

The following types of luminosity profiles are used for the lu- 
minosity decomposition: 



3.1.1. Disk 

For the disk, an exponential profile with the surface density 
distribution 



S_ED — So • e '■d 



(6) 



is used with Sg as the countrate in the center and the scale 
length rd determining the slope of the luminosity profile. Due 
to its inclination the disk appears elliptically projected with an 
apparent axis ratio b/ad and a position angle PAd of its appar- 
ent major axis. Free parameters therefore are: 

- the central intensity Eq 

- the scale length r^ 

- the apparent axis ratio b/ad 

- the position angle PAd 

The apparent axis ratio b/ad (ie., the inclination i) defines - 
together with PAd - the LON for the disk that will be used as 
deprojection axis for all three components. 



3.1.2. Bulge 

The central part of NGC 2336 is modelled by a profile of gener- 



alized Hernquist type (see Dehnen 1993) with the density pro- 
file 



Pq{m) 



with 



(3 - 7)M n 



1 



Anq ml (m + rf,)4-7 



2 ^ 

y + — 



(7) 



(8) 



1 results in a spherical 



q ~ yields the flat disk case, q 
distribution. 

With rb as a scale length and M as the integrated mass, the 
central slopes of the densities are oc r~^. This family of pro- 



files include several special cases: The steep Jaffe (1993) type 
(7 = 2), and the Hernquist ( |1990| ) type (7 = 1.0). Moreover, 
this profile type is a good approximation for the classical de 
Vaucouleurs (1948) profile (7 = 1.5). Free fit parameters are: 



- the scale length rb 

- the total mass M 

- the profile parameter 7 

- the vertical flattening of the bulge q 

- the apparent axis ratio b/ab 

For the position angle of the bulge, PAb, the same value as for 
the disk is used, since in disk galaxies the rotation axes of bulge 
and disk generally coincide. Nevertheless, b/ad and b/ab may 
differ due to the different thickness of the components. 

The projected luminosity density of the bulge that can be 
observed in the sky is obtained by evaluating the integral 



S(r) - 2 r 
Jr 



p{r)r 



yr 



i?2 



dr, 



(9) 



and is expressible in ter ms of e lementary functions for integer 
7 only (see Dehnen et al. 1993 for further details). 
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3.1.3. Bar 



A Ferrers (1877) bar with the surface density profile 



S(m^) 



So(l-m2) 



2^n 







ifm^ < 1 
ifrn^ > 1 



with 

m = 



(- 



a>b> 



(10) 



(11) 



is used. The total mass of the 2D bar obtained by integration is 

Mf,ar = T^ab^o^^r^ (12) 

r(n + 2) 

(cf. Matthias 1993). Free parameters that have to be adjusted 
during the fitting procedure are: 

- the apparent bar major axis a 

- the apparent axis ratio b/a 

- the central surface density Eq 

- the exponent n that determines the density distribution of 
the bar 

- the position angle of the bar major axis, PA^ar 

In our models, all Ferrers bars will be used with n — 2 for 
simplicity reasons. 

3.2. Results 

The model parameters are adjusted iteratively in order to min- 
imize the residuals that remain after the subtraction of the 
model from the observed luminosity distribution. The model of 
NGC 2336 that fits the observed light distribution best is shown 
in Fig. pi together with the remaining residuals. The parameter 
values obtained from the fit are given in Table 0. 

As can be seen from the residuals in Fig. H modelling the 
central regions is difficult: residuals in the bulge area are still 
clearly visible after the subtraction of the model. This is a con- 
squence of the fairly poor spatial sampling of the bulge area in 
the NIR images. In the same figure we notice that the spiral arm 
pattern is still present, because an average disk was subtracted 
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n = 1.58 


a = 4.17 


scale length [arcsec] 


42.9 


14.2 


37.5 


apparent h/a 


0.63 
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175° 
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120° 


relative contribution 


66.0% 
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profile parameter 7 




0.8 





Table 3. Parameter values for the decomposition model of 
NGC 2336. PA denotes the orientation of the major axes of 
the three components (counter-clockwise). Due to its intrinsi- 
cally rounder form compared to the disk, the bulge exhibits a 
larger b/a ratio although the major axes of both components 
are oriented in the same direction. 
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Fig. 3. Image slice through the three central pixel columns (up- 
per plot) an rows (lower plot) of the original J -band-image 
(solid line), the model (dotted line) and the remaining residu- 
als (dash-dotted line). For further explanations see text. 

generating positive residuals in the spiral arms and negative 
ones in the interarm regions. 

A quantitative impression of the quality of the fit is pre- 
sented in Fig. pi where the three central pixel columns and rows 
of every image in Fig. are averaged. The overall fit to the lu- 
minosity distribution of NGC 2336, though revealing small dif- 
ferences between the model and the observations, is sufficiently 
good. 

The morphological decomposition yields the result that 
NGC 2336 is dominated by the disk component, which con- 
tributes « 66% to the total luminosity, with the remaining light 
coming roughly equally from the bulge and the bar component. 

3.3. Deprojection 

Because of the incUned position of NGC 2336 in the sky, the 
2D model of the surface brightness has to be deprojected in or- 
der to obtain the intrinsic bar axis value, b/ai,ar,intr- We use 
the deprojection angles 0, 6, and ip, according to the definition 
given by Goldstein( |1981| ), with 6 being equivalent to the incli- 
nation i, and ^ representing the offset between at, 



'ar, app 



the LON. (p is only necessary if a non-axisymmetric compo- 
nent (a bar) is present and is used to give the bar the desired 
orientation before deprojectingthe disk. The results of the de- 
projection are shown in Table 0. The large difference between 
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Table 4. Deprojection parameters for NGC 2336. The Euler 
angle 9 is equal to the observed inclination of the galaxy, i, 
while the Eulerian ip denotes the bar offset from the LON. 



the apparent and the intrinsic bar axis ratio in Table Q originates 
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from the orientation of the bar in NGC 2336 with respect to the 
LON: Qbar is nearly orthogonal to the LON around which the 
galaxy plane is rotated during the deprojection, so the minor 
axis of the bar remains unchanged. The enlargement of atar 
then leads to a decrease of b/abar, intr- 

Since the deprojection angles are three additional free pa- 
rameters, the models will have to be tested for possible changes 
in the predicted kinematics that are caused by errors of (p, 9, 
and i/j. The most sensitive parameter is probably tp: It does not 
only include a possible misalignment of the LON, but also a 
wrong determination of the bar offset from this line. In case of 
NGC 2336, the latter is important due to the orientation change 
of the bar isophotes towards larger radii (isophotal twist). To 
keep the number of free model parameters as small as possible, 
(j) and 9 will remain fixed for all models with their photometri- 
cally derived values being used, tp will be varied independently, 
the results are presented in Section 10, where the effects of pa- 
rameter changes are discussed in detail. 

A last correction has to be made for the projection effects 
of the flattened spherical bulge of NGC 2336. An infinitesi- 
mally thin disk that is viewed pole-on represents an intensity 
/ that changes to Ig = 1/ cos 9 when viewed under the incli- 
nation 9. In contrast, a spherical bulge would look the same 
from every direction. It is therefore possible that a moderately 
flattened bulge introduces further uncertainties in the conver- 
sion factor Cbuige defined in Section 3 which converts coun- 
trates into luminosities. These uncertainties would be caused 
by changed inclination corrections for a non-spherical bulge. It 
will be discussed in section 10 that the exact shape of the bulge 
does not affect the kinematics of the NGC 2336-models signif- 
icantly, since for reasonable flattening ratios q = 0.4 ± 0.2 the 
projected mass density in the disk plane exhibits only minor 
changes. 

4. Kinematical Observations 

Spectroscopic observations were performed to determine the 
velocity fields of stars and gas in NGC 2336. By means of 
longslit spectroscopy of high spatial and moderate spectral res- 
olution the radial velocities of stars and gas clouds were mea- 
sured up to the outer regions of the disk (r > 100") with high 
S'/A^-ratios. 

Emission-line kinematics: 

For the kinematics of the warm HII gas component we used the 
light from the [Nil]- and i/a -emission lines in the red wave- 
length range. Since NGC 2336 is an intermediate-type galaxy, 
the disk contains enough gas for emission-line-spectroscopy 
even in the outer parts. Radial velocities were obtained by 
cross-correlating single rows of the longslit spectrum with the 
extracted central line from the core region of NGC 2336. The 
velocity dispersion of the warm gas lies far below the instru- 
mental resolution and was therefore not studied. 
Absorption-line kinematics: 

The spatial distribution of the stellar velocities and velocity dis- 
persions was examined by using absorption line spectra in the 
green wavelength range (4500 — 5500A). Using strong absorp- 



tion lines Hke Mgl at AArest=5167A, 5173A, 5184A, Cal fines 
at AArest=5262A, 5270A, and several Fe lines, the longslit 
spectra of NGC 2336 were cross-correlated with the spectrum 
of suitable template stars (e.g. of type G8III, KOI II, Kill I, 
K5III). The stellar velocity dispersions are much higher than 
those of the warm gas, they exceed the spectral resolution of 
our instrument configuration and can be studied in greater de- 
tail, at least in the central regions. 

All spectra were obtained during an observing run in March 
1997 using the MPIA 3. 5m- telescope at Calar Alto (Spain) 
with the TWIN-spectrograph. Due to a beamsplitting mirror 
with a crossover wavelength of Acr = 5500A, this device 
was capable of observing the blue and the red channel simul- 
taneously. The cameras in both channels were equipped with 
2048 X 1024 SITE-CCDs with a pixel size of 15Aim. With a 
scale length of 1" = 178/im along the slit and the correction 
for the ratio fcoii/ fcam — 6.34, we obtained a spatial reso- 
lution of 0.56"/pixel. Gratings with a dispersion of 36A/mm 
were used, the grating angles were adjusted such that the cen- 
tral wavelengths became Ac. blue = 5100A in the blue channel 
and Ac, red — 66OOA in the red channel. 

To achieve a S/N ratio as large as possible, the slit width 
was adjusted to 3.4". With our instrument setup and grating 
angles, the corresponding slit image was 4.5 pixels in the blue 
channel and 4.05 pixels in the red one. 

NGC 2336 was observed in 4 slit positions along the major 
and minor axes of disk and bai- {adisk, bdisk, abar, bbar)- The 
corresponding position angles were Pv4 — 5°, 94°, 28°, 118°. 
Exposure times were the same for emission and absorption 
spectroscopy: 2 x 3600sec in each slit position. The data re- 
duction process for both, the emission and the absorption line 
spectra, used a standard pre-reduction procedure containing the 
subtraction of the bias, the time-dependent linear dark current, 
and a flatfield division. The correction for the light distribution 
along the slit (including vignetting effects of the camera optics) 
was made by means of a slit profile obtained from a series of 
skyflat exposures. 

The subsequent main reduction process of the stellar spec- 
tra consisted of the following steps: Removal of saturated pix- 
els (cosmics), wavelength calibration (equidistant in Vr oc 
AA/A) via a rebin in log A, sky subtraction, continuum sub- 
traction, etc. . In a last step, these spectra were compared with 
correspondingly prepared spect ra of the template stars. We 
used a hybrid method of Bender (1990) which combines cross- 
correlation and Fourier quotient evaluation . This method is less 
sensitive to template mismatching and does not necessarily as- 
sume Gaussian broadening functions. From this procedure, the 
stellar radial velocity curves were obtained. Velocity disper- 
sions were computed using the broadening of the maximum of 
the correlation function. 

The data reduction of the emission line spectra was roughly 
the same, except for the fact that no template star, but the cen- 
tral line of the galactic spectrum was used for correlation, and 
for the slightly different correlation program being used. 

To collect as much light as possible, we chose a slit 
width of 3.4" for all observations which of course affects 



K. Wilke et al.: Mass distribution and kinematics of tiie barred galaxy NGC 2336 



the spectral resolution: With 0.54A/pixel (blue channel), the 
nominal resolution for measuring velocity dispersions was 
a = 31.7km/sec/pixel. Projected onto 4.5 pixels (blue chan- 
nel), the slit width became cr = 4.5pixelx31.7km/sec/pixel 
« 140km/sec in velocity space. But since we obtained the ve- 
locity dispersions from measuring the broadening of the max- 
imum of the correlation function only, the achievable resolu- 
tions are much better: A lower threshold of aeff ~ 65km/sec 
obtained for our stellar spectra is only marginally worse than 
the Nyquist-limited spectral resolution of slit. The resulting ro- 
tation curves of gas and stars and the stellar velocity dispersion 
curves are shown in Fig. Q, q| and g. 




r [arcsec] 

Fig. 4. Gas rotation curves of NGC 2336 with a spatial resolu- 
tion of 0.56" /pixel. For each of the 4 slit orientations the po- 
sition angles (PA) are given. The corresponding slit positions 
in the sky are indicated by the labels below the rotation curves 
and by the thick line in the isophotal plots (J-image). Due to 
the limited slit length there is a cutoff e.g. at PA — 5° at radii 
> +105". Notice that the S/N-ratio is remarkably high, even 
in the outer regions of the disk (at 100" j. 



Considering the gas rotation curve for a slit orientation near 
the kinematical LON (Fig. H, PA — 5°), we notice the steep 
rise of the radial velocities to values of \vrot\ ~ 185km/sec 
within 13" only, corresponding to 1.45kpc. This ascent in the 
inner part of NGC 2336 is caused by a centrally condensed 
and moderately massive bulge component. The presence of the 
bar is revealed by small humps (overshootings) in the rotation 
curve at a distance of r sa 13" from the center, after which the 
velocities decline again (|Atirot| ~ — 45km/sec). At distances 
from the center greater than 33", Vrot is constant, indicating 
that NGC 2336 does not contain a disk component with a large 
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Fig. 5. Stellar rotation curves of NGC 2336. Legends are the 
same as in Fig. ^ The S / N -ratios of the stellar rotation curves 
are lower than those of the gas rotation curves. The general 
kinematic behaviour of the stars differs significantly from that 
of the gas (for further details see text). 



scale length that would cause a considerable increase of Vrot 
beyond the regions dominated by the bar. 

In contrast to that, the stellar rotation curves for the same 
slit orientation {PA — 5°) in Fig. @ show a different behaviour: 
The final values of Vrot are not reached within of r « 24" 
(2.7kpc), although the final values of Vrot are indeed the same. 
Therefore the central slopes of these two gas rotation curves in 
the innermost regions are significantly shallower than those of 
the stellar rotation curves. 

This behaviour originates from the quantitatively different 
orbits which stars and Hll-clouds are moving on: Because of 
the non-neglectible cross-section of the spatially extended HII- 
regions, they accumulate preferentially on non-intersecting or- 
bits in the disk. If Hll-regions would populate self-intersecting 
orbits with loops, collisions would remove them from these 
paths after only a few rotation periods of the bar As a result, 
intersecting orbits do not play any role for the kinematics of the 
warm gas component. In other words, the radial component of 
the motions of the Hll-regions will be dissipated with time so 
that only circular or nearly circular orbits survive. This means 
that - viewed from the corotating reference frame of the bar 
- in our stationary models the Hll-clouds will populate closed 
non-intersecting orbits of elliptical or circular shape. 

In contrast to that, stars can be regarded to as point-like test 
particles that do not collide on timescales of H^^. In general, 
stars will therefore also populate orbits with a high contribution 
of radial motion. This is especially true for the stellar bulge 
which - like an elliptical galaxy - mainly consists of disordered 
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Fig. 6. Stellar velocity dispersions of NGC 2336. Legends are 
the same as in Fig. ^ Only the velocity dispersions of the bulge 
reach values well above the instrumental resolution. The disk 
cannot be detected. 



motion. Therefore the integration over many stellar orbits along 
the line of sight results in a reduced radial velocity and high 
velocity dispersion values in the bulge region (cf. Fig. ^. This 
causes the different slopes of the stellar and Hll-rotation curves 
in the inner regions of NGC 2336 in Fig. Q and g. 

An even larger difference between the motion of stars and 
gas occurs in case of a slit position along the minor axis of the 
disk {PA = 94°) in Fig. ^ and g. The humps and dips in the 
gas rotation curves with jAwrotl = llOkm/sec at a distance of 
r = 7" (0.6kpc) from the center indicate high streaming veloc- 
ities along the bar Those humps and dips are also present in 
the stellar rotation curve in Fig. g, which is a necessary condi- 
tion for the existence of the bar structure, since stars constitute 
the bar But as is clearly visible, those overshootings are much 
smaller, since the innermost regions of NGC 2336 are domi- 
nated by the highly disordered radial motions of the bulge stars 
which nearly cancel out during the integration along the line of 
sight. To summarize, Hll-clouds can been regarded to as ideal 
test particles to trace the underlying gravitational potential of 
the bar The latter is mainly generated by the stars, carrying 
most of of the total visible mass in normal intermediate-type 
disk galaxies. 

The comparison between radial velocities of gas and stars 
for PA = 118° (along the major axis of the bar) in Fig. 
and p] yields no qualitative difference, except for the 5/A^-ratio 
being again low for the stellar kinematics, which prevents us 
from comparing the outer regions with r > 50". Both curves 
(gas and stars) show similar solid body rotation over the whole 
length of the bar 



As shown in Fig. pi NGC 2336 reaches a central velocity 
dispersion of ct = 160km/sec. In all velocity dispersion curves 
the bulge extends up to radii of r « 40". This means that the 
apparent bulge axis ratio cannot be much smaller than 1, with 
b/a « 0.8 as a lower boundary, taking the velocity dispersion 
curves at PA = 5° (disk major axis) and PA = 94° (disk 
minor axis) as reference. It is not possible to determine this axis 
ratio here more exactly, but a kinematically derived condition 
b/a > 0.8 agrees well with the photometrically derived value 
of6/a = 0.85. 

5. Potentials 

To derive the common potential of disk, bulge, and bar, we 
make the following approximation: Disk and bar are collapsed 
to their plane of symmetry, their potential is calculated in 2D. 
The bulge is considered as an oblate spheroid, but the poten- 
tial of this 3D mass distribution is computed only in the disk 
plane. The following section deals with the derivation of the 2D 
potentials from the morphologically derived luminosity profile 
and their implementation into the FORTRAN-/C-orbit integra- 
tor. 

In general, potentials and density distributions are con- 
nected by Poisson's equation A$(a;) = AttGp{x). For the 
Ferrers bar, the exponential disk and the generalized Hernquist 
bulge, no analytical closed expressions for their 2D potentials 
are available (except for some special cases), so they have to be 
evaluated by means of numerical representations or expansion 
series. 

5.1. Disk 

It is not appropriate for the potential calculation to consider 
the disk as a collapsed 3D sphere with thickness zero, since 
one encounters double integrals that cannot be solved easily. 
As a better approach, we use Bessel functions t o deri ve the 
disk potential (described by Binney and Tremaine 1987| , origi- 
nally developed by Toomre 1962 ). With the ansatz $fe(-R, z) = 
e~'"'l^l Jo(fci?), the Laplace equation V^<I> = can be solved 
except for the 2D plane of the disk, where the gradient is 
discontinuous. Using the Gauss theorem, the mass density 
that generates the gradient discontinuity is found to Sfc = 
— 2^ Jo(fc-R) with Jo as cylindrical Bessel function of 0th or- 
der With the Hankel transform S{k) defined by 



(13) 



I](fc) = / S{k)EkiR)dk, 
Jo 

the potential is given by 

/■OO 

$(i?,z) = / S{k)^k{R,z)dk. 
Jo 

It can be shown that for the exponential disk with 

S(i?) = Eoe-^/^^ S{k) takes the form 



(14) 
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For the 2D disk (z — 0) and with y = R/2Rd, this results in 
$(i?,0) - ~nG^oR[h{y)Ki{y) - h{y)Koiy)]. /„ and K^ 
are modified Bessel functions of 1st and 2nd kind. This expres- 
sion is implemented in the program used for the calculation of 
closed particle orbits. 

5.2. Bulge 

For the generalized Hernquist bulge luminosity profile, there 
are no closed analytic expressions for the corresponding $(r) 
(see Dehnen 1993). Therefore, our orbit integrator includes 
FORTRAN-routines written by Walter Dehnen (Oxford) which 
allow a numerical treatment. The 3D bulge potential for the 
generalized Hernquist profile takes the form 
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The bulge is the only component with a 3D potential imple- 
mented into the numerical orbit integrator programme (disk 
and bar are considered with the corresponding 2D expressions). 
Since the flattening parameter q determines the mass density 
projected onto the 2D plane, ^buige will be sensitive to changes 
of that parameter. 

5.3. Bar 

To illustrate the derivation of the bar potential ^bar, we first 
derive a 3D expression for the potential of a general oblate 
spheroid and reduce that expression to the 2D case afterwards. 
The resulting 2D potential allows a very simple numerical 
treatment of the bar. 

A generalized 3D-disk can be expressed by means of an 
expansion method described by Binney and Tremaine ( 1987 ) 
that uses oblate spheroidal coordinates. Those coordinates use 
the (j) from cylindrical coordinates, but replace R and z by 
R — E ■ cosh u sin v, and z — E ■ sinh u cos v, respectively, 
with S as a constant parameter determining the range of the 
coordinate system. The potential of a 3D disk outside the inner 
2D area u = 0, where this coordinate system would be am- 
biguous (there exist two w-values for each R), takes the form 
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Q™ are the Legendre polynomials of 2nd order, Y" 
spherical harmonics. Vim is a constant that contains 
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as the generating surface densities that account for the ^im- 
gkrn and G are constants. The generating 3D potentials 4>/m 
are divergent for i? ^ S, but S remains finite. This causes 
numerical problems when one wants to compute potential val- 
ues at i? = S. This restriction limits the use of the expansion 
series approximation, because we are prevented from calculat- 
ing <1> at large radii (in the vicinity of S). Nevertheless, the 
method is generally useful for approximating potentials with 
sharp boundaries R{ip) (like the Ferrers bar). 

In the inner 2D disk (defined by u = in the spheroidal co- 
ordinate system), which is the only interesting region for calcu- 
lating the potential of a 2D bar, the coordinate transformations 
reduce to i? = S • sin v and z — Q with < u < 7r/2. The 
collapsed 2D potentials are no longer u-dependent: 



*im = ViynYim{v,Lp). 



(23) 



This simple form is sufficient for our 2D studies and has the 
advantage that arbitrary bar forms (e.g. asymmetric or drop- 
shaped bars) can be treated. 

According to the initial definition as scale parameter, S de- 
fines the model unit for the radius, so S = 1 is the maximum 
radius at which the potentials $;„ can be calculated. Scaling of 
the model according to different bar lengths or different radii 
up to which the models have to be examined is achieved by 
defining abar and hbar as fractions of S = 1, e.g. a = 0.7 and 
h = 0.3. Scale lengths of the bulge and the disk component 
are then adjusted according to the bar length. Therefore, in a 
model with a bar major axis a, potentials can be calculated up 
to a distance of 1/a. For all models a bar length abar = 0.7 will 
be used, the scalelengths of disk and bulge are then defined ac- 
cording to abar ■ It turns out that, with the resulting maximum 
radius of 1/a = 1.42 bar lengths accessible by the potential 
expansion series, the outer Lindblad resonance OLR lies well 
within the reach of most of our models. To restrict the con- 
sumption of computational time, we truncate the expansion se- 
ries (equation 20) at appropriate values m,,nax and Imax that 
determine the angular and radial resolution. For most of the 
models discussed below we use Imax — 50 and rumax — 12 
which ensure sufficiently resolved potentials. 

In a final step, the total potential of the NGC2336-models 
is obtained by simply coadding all contributions: 



(20) $ 



total — ^disk + ^bulge + ^bar- 



(24) 
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6. Resonances and periodic orbit families 

For small deviations from axisymmetry, barred potentials can 
be treated by means of linear theory, in which a star radially 
oscillates around its guiding center with an epicyclic frequency 
K. The guiding center moves around the center of the galaxy 
with the circular frequency n{R). In this approximation, sev- 
eral resonances appear where particles are subject to strong in- 
teractions with the rotating bar potential: the inner Lindblad 
Resonance (ILK), the Corotation (CR), and the outer Lind- 
blad Resonance (hereafter OLR, for an extended review see 



Sellwood and Wilkinson 1992 ). However, in case of NGC 2336 



this approximation is not valid due to the non-neglectible con- 
tribution of the massive bar component to the total potential. 
But although single particle orbits can no longer be computed 
analytically, the examination of the resonances is still useful 
when studying the kinematics. The locations of the resonances 
have to be determined directly from the analysis of the phase 
space (Combes 1994). 

In all NGC2336-models, the typical xi, X2- and X4- 
families of closed, therefore periodic orbits exist between the 
different resonance locii. The general properties of these or- 
bit families are the same for all NGC 2336-models and repre- 
sentative for a disk galaxy with a bar component contributing 
roughly 10% to the total potential, as given by Contopoulos 
( [1980| ). 

The orbits of the xi-family are elongated along the bar 
major axis within CR, they are the main orbit family sup- 
porting the bar. They cover the range from the center up to 
the outer regions of NGC 2336. Outside CR they become 
rounder and change their rotation sense from prograde to ret- 
rograde, since particles on such orbits are moving slower than 
the bar. Epicyclic orbits with loops, belonging to higher (n:m)- 
resonances of the xi -family, will be populated by stars only 
and are therefore excluded from the orbit sets used for the 
construction of artificial rotation curves of the Hll-gas. 2:2- 
orbits are confined to the region between the center and ILR. 
The nearly circular a;4-orbits do not support the bar. Because 
of their retrograde orientation they are not populated by gas 
clouds. They are therefore neglected for the construction of 
Hll-rotation curves. 

In the rotating reference frame of the bar, the normal en- 
ergy E is no longer conserved along the particle trajectories, 
since virtual forces appear that add a centrifugal and a Corio- 
lis term to the equation of motion. Instead of E = ir|„ + $ 
(with Tin denoting the velocity in an inertial coordinate system, 
cf. BT87), the so-called Jacobi energy Ej is conserved: 



dEj 
dt 



-0, with Ej^-r^ + ^--\npXr\\ 



(25) 



Ej will therefore be used as model parameter instead of the 
normal energy E. 



7. Model Units 



7.1. Length Scale 

As described in Section 5.3, model length scales are coupled to 
the numerical expansion radius S = 1 by choosing the maxi- 
mum radius the kinematics of orbits shall be computed for. It 
turns out that for reasonable mass distributions, 0.7S is a good 
choice for abar, since the OLR lies then well within the reach 
of our models. The length units LU are calibrated using the 
NIR-images. As given in Section 3.2, the best photometric 
decomposition model of NGC 2336 has a minor axis length of 
the bar of b ^ 27.0 pixels. With the scale of 0.67'7pixel of the 
MAGIC camera, we get b = 18.09". With a b/abar, mtr = 0.28 
the length of a becomes 

a ^ 18.09'70.28 = 64.6" = 0.75. 

With a distance of d = 22.9Mpc (based on Hq = 
lOOkm/sec/Mpc), we get 1" = lllpc, which results in the 
model scale (LU) 



ILU = 92.3" ^ 10.25kpc, 



(26) 



so the photometrically derived bar of NGC 2336 has the dimen- 
sions a ~ 7.2kpc and b ~ 2.0kpc. 

7.2. Mass Scale 

We calibrate the mass scale via the equation for the circular 
velocity 



d(j} _ CM(r) 



(27) 



"^ dr r 

We use those velocity values that are found far out in the disk 
(r K, 9kpc) to obtain 



Mtot = Mbulge + Mdisk + Mbar 
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(28) 
(29) 



Since the relative mass-to-light ratios Cdisk and Cbar have to 
be found by optimizing the kinematical model, this is an itera- 
tive process. The kinematical effects of varying Cdisk and Cbar 
are studied by corresponding model sequences, similar to those 
for the scale lengths, the pattern speed of the bar, etc. With this 
method, the mass unit MU is determined to 



1MU= 3.234- IO^Mq. 



(30) 



7.3. Time scale 

The time scale is adjusted such that G = 1. This yields 

ITt/ EE 2.71 • 10 Vs- (31) 

8. Phase Space Analysis and the Integration of Closed 
Orbits 

8.L The Poincam-Surfaces of Section 

Stars and gas particles in real barred galaxies are normally 
not likely to strictly follow periodic orbits (so their occupation 
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number should be low), but it can be shown that most of the 
non-periodic orbits in barred galaxies are trapped to oscillate 
about a parent periodic one (for an extended review on phase 



space analysis see Sellwood and Wilkinson 1992). Such orbits 
are often referred to as quasi-periodic. 

Therefore the examination of periodic orbits is an impor- 
tant tool, because the structure of the bar is largely shaped by 
those parent orbits. To derive the spatial extent of the regions 
supported by the parent orbits, the division of phase space has 
to be studied by means of numerical methods. For those pur- 
poses, cuts through phase space at a certain orbit energy - 
the Poincare-Surfaces of Section (hereafter SOS) - are widely 
used. Positions and velocities are measured every time a parti- 
cle crosses the plane x — 0, x > (or y = 0, y > 0), thereby 
constituting the y, y- (or (x, i))-Surface of Section at a given 
Jacobian energy Ej. 

The division of phase space in the case of the NGC 2336- 
models exhibits characteristic features which are explained us- 
ing results of a model with a moderately flattened (q = 0.6) 
bulge (model parameters are listed in Table |]). The (x, x)- 



bar 


atar [LE/kpc] 
btar [LE/kpc] 
b/a 

^bar 


0.7/7.17 
0.21/2.15 

0.3 

1.5 


bulge 


7 

rt [LE/kpc] 

g 

Obulge 


0.9 

0.20/2.04 

0.6 

1.0 


disk 


ra [LE/kpc] 

Cdisk 


0.46/4.76 

2.5 


Qp [km/sec/kpc] 


17.85 



Table 5. Parameters for a NGC 2336-model with a moder- 
ately flattened bulge and low C-values for bar and disk. All 
scale lengths are given in model (LU) and physical units (kpc). 
Masses are given in solar mass units (Mq). The mass calibra- 
tion refers to a preliminary calibration according to Sect. 3. 
With the selected flp-value, the CR is placed at the end of 
the bar: rcR ~ 1.0af,ar = 0.7LU. With a bar axis ratio of 
b/a — 0.3, a value slightly larger than that obtained from de- 
projection (b/a = 0.28) is chosen. Ej is computed according 
to its definition in Section 6, using the model units of Section 7. 
The same units are used for all following tables. 



Surfaces of Sections for different Ej-\alues are displayed in 
Fig. ^ For low _Ej-values (left plot in Fig. ^, the phase space 
of the NGC 2336 is dominated by xi- and a;2-orbits on the 
prograde side (x > 0) and by a;4-orbits in the retrograde half 
(x < 0) of the Surface of Section. Each periodic orbit is sur- 
rounded by closed invariant curves (hereafter ICs) of quasi- 



periodic orbits. For this energy well below that of the inner 
Lindblad resonance {EiLp/), a::2-orbits cover a considerable 
area in phase space. This behaviour originates from the cen- 
trally condensed model with a high bar mass that produces a 
strong ILR. Notice that particles at that i?j-values can travel a 
small region within the bar only, since the S'OS'-boundary ex- 
tend only up to |a;| w 0.05LC/, while the bar axes lengths are 
a = 0.7LC/ and b = 0.245LC/, respectively. 

At higher i?j-values (central plot in Fig. 0), the retrograde 
side remains nearly unchanged, while on the prograde side 
a;2-orbits are now completely absent since Ej > Ejlh. The 
xi -orbits have become rounder and a small amount of semi- 
ergodic motion can be traced that surrounds the a;4-island. 
These semi-ergodic orbits, indicating irregular motions, are al- 
ways present in models with a massive bar and hint at a self- 
regulating mechanism for the bar strength: The more mass is 
accreted in the center of a barred galaxy, the steeper the cen- 
tral potential will become and the higher will be the percentage 
of semi-ergodic orbits. But with too much semi-ergodic mo- 
tion, a bar will not be able to sustain its length and mass on 
longterm scales because of a lack of particles on orbits support- 
ing the bar. If therefore two NGC 2336-models with the same 
artificial Hll-rotation curves contain different amounts of semi- 
ergodicity at comparable energies, the model with the smaller 
amount of semi-ergodic motion will be chosen for stability rea- 
sons. 

Slightly above the CR energy Ecr. = —26.5, the Surfaces 
of Section look completely different. At Ej = —26.3, the re- 
gions outside the bar become energetically accessible to closed 
orbits, as can be seen in Fig. ^ (right plot). But the phase space 
outside CR is only partially populated with orbits, which is 
a tendency that becomes even more evident when proceeding 
towards even higher energies. For this model, CR is placed 
at TCR — 0.7 (at the end of the bar) by adjusting the pattern 
speed Up. In the prograde half of the phase space, numerous 
islands belonging to (n:m)-families of higher resonant orbits 
have replaced the xi-orbit and its surrounding ICs. For the 
Hll-kinematics, those higher resonances are of no importance 
due to their small occupation number, which is caused by their 
self-intersections leading to a depletion after short times. They 
contribute only to the stellar bar in a significant manner. The 
most important of those islands belongs to the (4: l)-resonance, 
which is located at a; w 0.28 in Fig. ^ (right). Other higher 
resonant families (e.g. the (6:1) family) are less important. 

8.2. Construction of Artificial Rotation Curves 

Our orbit integrator uses Hll-clouds as test particles to trace 
the underlying gravitational potential ^totai of the NGC 2336- 
models. Trajectories are integrated on a grid-based scheme, 
which means that in each cell of the phase space Hll-clouds 
with certain initial conditions (e.g. energy, velocities) are 
started and their trajectories are followed through the phase 
space. 

For each model, complete sets of closed orbits that belong 
to the xi- and (if they are present and/or necessary) a:;2-family 
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Fig. 7. (x, x)-Surfaces of Section for a typical NGC 2336-model with a moderately flattened (q ~ 0.6j bulge. The Surfaces of 
Section are computed for Jacobi energies of Ej = —45.0 (left), Ej = —30.0 (center) and Ej — —26.3 (right). Below every 
SOS, the periodic orbits of the main families being present in the SOSs are plotted in the (x, y)-plane, the solid lines indicating 
the positions those orbits can be found at in the SOSs above. 



are computed: The orbit integrator uses the phase space coor- 
dinates of the periodic orbits at various i<^j-values as starting 
points and computes a number of orbits (normally 100 of each 
family type) belonging to different particle energies. The pro- 
gram uses a continuous check for the particle energy during 
the integration of every orbit to ensure conservation of the Ja- 
cobian energy Ej . The energy interval for which the orbit sets 
are computed is adjusted with regard to the different Lindblad 
Resonances in the model (e.g. orbits between ILR and CR or 
CR and OLR). 

All resulting orbit sets are then projected to the sky with 
the same inclination and orientation as the galaxy. The orbits 
generate the kinematical "backbone" of our artificial galaxy. A 
virtual longslit is projected to the sky as well. Corresponding 
to the observations, a constant slit width of 3" is chosen for 
all models discussed below. The orbit integrator looks for the 
points where an orbit crosses the virtual slit and computes r and 
Vrad for cvcry crossing point. All crossing points together yield 
the artificial rotation curves for a given orientation angle of the 
virtual slit. The large number (> 100) of orbits in each com- 
plete set ensures a sufficient spatial resolution of the artificial 
rotation curves even in the central regions of the NGC 2336- 
model. 

In a last step the rotation curves are convolved with a gaus- 
sian point spread function of constant width (1.5") to correct 
for the seeing conditions of the original longslit spectra. By 
varying the position angle of the virtual slit, all position angles 
from the observations can be traced. 



9. Variations of the Model Parameters 

The model parameters were varied to check for qualitative and 
quantitative changes in the synthetic kinematics of the HII gas. 
We also tried to answer the question whether a possible best- 
fitting model could be subject to ambiguities that are caused by 
the fact that the variation of two independent model parame- 
ters may have the same effect on the artificial rotation curves. 
Therefore, the following parameters were varied in a range that 
agrees well with the photometrical observations: 

1. The bar offset from the line of nodes, i/i: Several cases 
were calculated with this parameter varied in the range 
if) — 102° . . . 106° which is well in agreement with the 
J-band-observations. Since the Euler angles and (j) ^s 
subject to very small deprojection uncertainties only, they 
remain fixed for all models. 

2. Major and minor axis of the bar, a and b: Due to the 
presence of another dominant component in the center of 
NGC 2336 - the bulge - errors may occur when evaluating 
b/a\jar- In principle, a larger h/abar could be compensated 
by a different choice of the bulge scale length, r^. 

3. The bar mass, M\jar'- The exact determination of M^ar is 
not easy since the bulge component contributes to the lu- 
minosity in the inner regions as well. The bar mass (equiv- 
alent: Char) was varied at the expense of the bulge and vice 
versa. 

4. The bulge scale length, r^: The fit of the bulge scale length, 
Tb, to the J-band luminosity distribution is difficult because 
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of the small detector area covered by this component. Vari- 
ations of Tb in the models are therefore necessary. 
The flattening of the bulge, q: This parameter is important 
because it determines the percentage of the total projected 
bulge mass that is concentrated in the disk plane. The value 
of q could in principle be estimated directly from the photo- 
metric data by evaluating the different apparent axis ratios 
of the (flattened) bulge and the (infinitesimally thin) disk, 
b/aapp, bulge and h/aapp.disk- But as it is the case for r^, 
the reliability of the estimated 6/aapp. buige-value suffers 
from the low resolution in the bulge area. In our models, q 
is therefore varied in the reasonable range q — 0.2 ... 0.6, 
well consistent with the photometrically derived apparent 
axis ratio. 

The bulge mass, Ivhuige'- The effects of different bulge 
masses (CfjuZge-values) strongly mix up with the results 
of changing M^ar (Char), therefore models have to be 
checked for this parameter indepenently. 
The disk scale length, r^: The exact determination of r^ 
is limited by the size of the iV/i?-detector that partially 
omits the outer disk regions. Since no further information 
about Trf can be obtained, the photometrically derived value 
is used for all models. 

The disk mass, Mdisk- Changing Mdisk (Cdisk) causes sig- 
nificant velocity changes in the regions far outside the bar 
Its contribution to the total mass varies according to the 
masses of bulge and bar. 



Model 


1 


2 


bar 1 


abar 
bbar 

b/a 

^bar 


0.7/7.17 

0.245/2.51 

0.35 

1.0 


0.87/7.17 
0.245/2.51 

0.35 

5.0 


bulge 1 


7 

Tb 
1 

^hulge 


0.9 

0.20/2.04 

0.4 

1.0 


0.9 

0.20/2.04 

0.4 

1.0 


disk 1 


rd 

Cdisk 


0.46/4.76 

2.5 


0.46/4.76 

2.5 


ftp 


18.2 


19.1 



Table 6. Model parameters for two models with different bar 
masses. The scale lengths are given in model (LU) and absolute 
units (kpc). All other parameters are the same as in Table m To 
illustrate the qualitative difference between low and high bar 
masses, the bar conversion factor is set to an (observationally 
not confirmed) value of Cbar = 5.0 in model 2, while model 
1 uses Char = 1, corresponding to the mass directly obtained 
from the morphological decomposition procedure. 



9. The bar pattern speed, ^p. This parameter does not only Section 9.2, the effects of two different disk masses Mdisk are 



depend on the radial mass distribution of a model, but can 
also be varied independently since no direct observational 
constraint is available. To make the models better compa- 
rable when varying the other parameters, each NGC 2336- 
model should have CR placed at the same absolute radius, 
since certain features in phase space are spatially connected 
to the CR radius, e.g. the regions of semi-ergodic motion 
or the occurrence of higher resonant families. To achieve a 
constant value for rcR, it is necessary to change ilp each 
time the radial mass distributions are changed by variations 
of scale lengths and masses. On the other hand, f2p has to 
be varied independently (with all other parameters fixed) to 
study the influence of different corotation radii rcR on the 
NGC 2336-kinematics. The results of those variations are 
helpful for positioning CR with regard to the bar length. 
This is a step that can only roughly be performed using 
morphological features (e.g. bifurcations, spiral arms or 
rings). 
10. The profile parameter of the bulge, 7: Different values of 7 
change the radial mass distribution of only the bulge com- 
ponent. The effects of such changes are already examined 
by the variations of r^ in a sufficient manner. Therefore, we 
adopt 7 = 0.9 for all models. 

To illustrate the effects of the variations of the model pa- 
rameters, the Surfaces of Section and the rotation curves of four 
different models are shown in the following two Sections. Sec- 
tion 9. 1 deals with two models with different bar mass Mbar- In 



studied. Section 9.3 describes the effects of various J7p-values. 

9.1. Variation of the bar mass Mbar 

The bar mass is changed by variations of Cbar- Table o lists 
the basic parameters for two models with different bar masses. 
Model 1 with a M / L of Cbar = 1.0 uses the observationally 
determined value. Model 2 uses Cbar — 5.0 which obviously 
is not confirmed observationally but is useful for illustration 
purposes. In both cases, CR is placed at the end of the bar: 
rcB, = 0.7Lt/. 

In general, the division of phase space resembles the one of 
the introductory model: /Cs around the periodic xi-orbit are 
dominating the prograde side and x^-ICs cover the retrograde 
one. a;2-orbits are completely absent for both models, even at 
very low _E,/-values that are not shown here. 

Comparing the Surfaces of Section in Fig. |[ a much larger 
percentage of semi-ergodic motion for model 2 (lower row) can 
be detected at every energy level. This result supports a well- 
known fact: massive bars with their strong non-axisymmetric 
potential favour the existence of numerous semi-ergodic orbits. 
In the SOS of model 2 at low energies (lower row, left plot), 
a small number of islands around the /Cs of the X4-family is 
located within a sea of irregular or semi-ergodic motions. Due 
to the limited resolution of the grid and the limited time over 
which the non-periodic orbits are traced by our integration pro- 
gram, the space between the /Cs is not completely filled wih 
islands and single points. 
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Fig. 8. SOSs for model 1 (Cbar — 1, upper row) and 2 (Ctar = 5.0, lower row). Energy values are increasing from left to right. 
The values are: Ej = —30.0, —24.8, —24.6 (model 1) and Ej = —40.0, —26.5, —26.3 (model 2). Corotation is placed at 
rcR = 0.7. 



Further differences between model 1 and 2 occur when 
one approaches higher energies (central and right columns 
in Fig. pb around Ecr = —24.7 (model 1) and Ecr = 
—26.4 (model 2), respectively. The (4:l)-resonance, visible at 
X ~ 0.25 . . . 0.3 (Fig. g[ central column), is much stronger for 
model 1 than for model 2. 

When comparing the artificial Hll-rotation curves of both 
NGC 2336-models in Fig. ^, significant differences occur 
for nearly all slit orientations: Model 2 (right column) pre- 
dicts higher velocities for the bar-dominated inner region of 
NGC 2336, especially for slit positions that trace the stream- 
ing velocities of particles which are moving on a; i -orbits along 
the bar (PA = 94°). While model 1 (left column) seems 
to be a good description for the Hll-kinematics for slit ori- 
entations along the major axis of the disk {PA = 5° and 
PA — 28°), model 2 produces velocities that are too high up 
to regions around r « 25". Only along the major axis of the 
bar {PA — 118°), no differences between model 1 and model 
2 are visible. 

Model 1 mainly suffers from the fact that the overshootings 
in the rotation curves at PA = 94° cannot be traced to their 
full extent: velocities are too low by | Aw^ot | ~ 45km/sec when 
comparing the peak values of observed and predicted rotation 



curves. In addition, the radial extent of these humps is too small 
by a factor of 2. 

A general mismatch between observations and model cal- 
culations occurs in the innermost regions up to r « 5", where 
the slopes of the artificial rotation curves of model 1 and 2 are 
too steep, compared with the observed ones. This is a hint at a 
too high bulge mass concentration in both models. 

The results of this section show that Cbar — 10 is too low 
to explain the observed velocity field in the inner regions, while 
Cbar — 5.0 is too high. The optimal value should lie in the 
range 1.0 < Cbar < 2.0. 

Second, 6/a-ratios which are close to the photometrically 
derived value (0.3) or agree with it, are reasonable, since it 
seems possible to achieve a good agreement between the model 
and the observations. If discrepancies between observed and 
predicted velocity field are still remaining for observationally 
confirmed bar parameters, other parameters which are closely 
connected to the kinematics of the inner regions of NGC 2336 
will have to be adjusted. E.g., the velocity contributions that 
are still missing in model 1 at PA = 94° with reasonable bar 
parameters could be generated by slight changes of fip or by 
raising the central mass {Mbuige)- 
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Fig. 9. Theoretical gas rotation curves of model 1 (Char = 1, 
left column) and model 2 (Char — 5.0, right column). Only xi- 
orbits are considered. Observed points are always marked by 
small dots without error bars. 



9.2. Variation of the disk mass Mdisk 

To examine the effects of different disk masses Mdisk, Models 
3 and 4 include disks with a Cdisk of 5.0 and 2.5, respectively. 
In contrast to models 1 and 2, CR is placed at rcR ~ \.lai,ar 
in both cases. 

The main difference between the two models is the pres- 
ence of a periodic a;2-orbit in model 4 at x « 0.09 (Fig. [lO| ), 
surrounded by /Cs. Those orbits are completely missing in the 
corresponding Surface of Section of model 3 at comparable 
£^j-values. Model 3 produces a::2-orbits only in the very cen- 
ter of NGC 2336. This fact is illustrated by the corresponding 
a;2-rotation curves shown in Fig. |l2| which cover a very small 
region only. Additionally, the closed /Cs of model 3 around the 
a;4-orbit are surrounded by islands and semi-ergodic regions al- 
ready at those low energies, a feature that is missing in model 
4. 

Semi-ergodicity becomes more important when we pro- 
ceed to higher energies (central and right column in Fig. nOb: 
At Ej = —32.4, in model 3 the prograde side of the phase 
space is dominated by a very strong (4: l)-resonance at a; ~ 0.2 
which is surrounded by large regions of semi-ergodic motions 
(upper row, central plot). In contrast to this, model 4 (lower 
row, central plot) exhibits a much weaker (4:l)-resonance at 



Model 


3 


4 


bar 1 


a-bar 
bbar 

b/a 

^bar 


0.7/7.17 

0.245/2.51 

0.35 

1.5 


0.87/7.17 

0.245/2.51 

0.35 

1.5 


bulge 1 


7 

Tb 
1 

Cybulge 


0.9 

0.20/2.04 

0.4 

0.63 


0.9 

0.20/2.04 

0.4 

0.63 


disk 1 


rd 

Cdisk 


0.46/4.76 
5.0 


0.46/4.76 

2.5 


ftp 


20.3 


22.2 



Table 7. Model parameters of models 3 and 4 with different 
disk masses. Scale lengths are given in model (LU) and physi- 
cal units (kpc), all other units are the same as in Table Notice 
the changes of Q.p that are caused by the variations of the ra- 
dial mass distribution. 



X w 0.3 and nearly no semi-ergodic motion at comparable 
energies {Ej = —35.8). Even at the highest energy values 
slightly above Ecb. (right column), semi-ergodic motion plays 
a less important role in case of model 4 than in model 3. 

Both models are not able to reproduce the observed kine- 
matics with sufficient accuracy: A Cdisk of 5.0 produces xi- 
velocities that are much too low for the central xi -orbits (model 
3, Fig. 1 1 left column). A moderately lower Cdi^fc-value of 2.5 



yields a better fit (model 4, Fig. ITTl right column, see also the 
results of the optimal model in Fig. 13 which were obtained 
with Cdisk — 2.5). In addition, model 3 suffers from the fact 
that the humps and dips in the observed Hll-rotation curves at 
PA = 5°, 28° and 94° are not reproduced to their full extent. 
This is caused by the large overall disk mass (Cdisk = 5.0), 
with a correspondingly small non-axisymmetric contribution of 
the bar to the total potential. Therefore, the xi-orbits support- 
ing the bar structure become less elongated, and the streaming 
velocities of the Hll-clouds along the major axis of the bar are 
much lower 

However, a simple reduction of Cdisk with all the other pa- 
rameters remaining fixed is no solution either, as is quite obvi- 
ous from the results of model 4 (Fig. O, right column): With 
Cdisk — 2.5, the central slopes of the rotation curves of this 
model are too steep due to the increased importance to the non- 
axisymmetric bar contribution to ^totai and to the small bulge 
scale length rf,. Additionally, the maximum velocities of model 
4 exceed the observed ones by |Aurot| ~ 50 — lOOkm/sec at 
PA = 5° and PA = 28°. Another feature of model 4 is the 
strong decline of the rotation curves after reaching the (correct) 
maximum velocities for a slit orientation of PA ~ 94°, which 
is not supported by observations, either. 
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Fig. 10. Surfaces of Section for model 3 (Cdisk — 5.0, upper row) and model 4 (Cdisk = 2.5, lower row). Ej-values are (from 
left to right): —40.0, —32.4, —32.2 (model 3) and —55.0, —35.8, —35.5 (model 4). In both models, Corotation is placed at 
rcR = 1-labar- 



It can be shown that a slightly larger bulge mass (Cbuige ~ 
1.0 instead of Cbuige = 0.64), increasing the mass concen- 
trated in the center of NGC 2336, extends the Hll-streaming 
motions along the bar to larger radii (PA — 94°). At the same 
time, a larger rt-value would induce a flatter rise of the veloci- 
ties in the innermost part of NGC 2336 and reduce the observed 
humps in model rotation curves at PA ~ 5° and 28° to the cor- 
rect height. 

Both models produce 2:2 -orbits, the corresponding rotation 
curves are also computed and shown in Fig. |l2[ While model 3 
(left column) generates 0:2 -orbits only in the very central region 
(< 3") of NGC 2336, in model 4 (right column) a::2-orbits are 
present up to radii of r w 15" in the case of PA = 5° and 28°. 
As far as those orbits are present in both models, no signif- 
icant differences occur between the a;2-rotation curves. How- 
ever, supporting the results of models 1 and 2, the a;2-orbits of 
models 3 and 4 cannot explain the kinematics of the inner re- 
gions of NGC 2336, which is clearly visible for a slit position 
of PA = 94°. To summarize, large Crf^sfe-values lead to an 
increased axisymmetric potential component, therefore the bar 
becomes less important. The consequences of the massive disk 
in model 3 are xi -streaming velocities that are in general too 
low by far to reproduce the observed velocity field. We draw 



the conclusion that a final model of NGC 2336 should use a 
Cdisk < 3. The necessary adjustment of the mass distribution 
in the central region could then be achieved by varying other 
parameters, e.g. the bulge mass, which was chosen especially 
low for the models 3 and 4. The slope of the rotation curves 
could be adjusted by varying r;,, the streaming velocities by 
lowering Up. 

9.3. Variation of the bar pattern speed fip 

In contrast to the parameters resulting from the morphologi- 
cal decomposition (e.g. component masses and scale lengths), 
Q,p cannot be observed directly. Therefore a sequence of iden- 
tical models with different fJp-values is examined. In general, 
all models produce similar velocity fields. Major differences 
occur only for a slit orientation of PA = 94°, where models 
with lower Q.p predict significantly larger humps in the rotation 
curves. This is a consequence of the fact that the rotation curves 
result from the xi -streaming velocities - which are identical for 
all models - subtracted by the various fJp-values. In contrast 
to that, the other slit orientations show only minor deviations. 
Obviously, variations of Vtp are an appropriate tool to adjust the 
velocities perpendicular to the bar major axis. 
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Fig. 11. xi-mtation curves of models 3 (Cdisk — 5.0, left col- 
umn) and model 4 (Cdisk — 2.5, right column). Observational 
values are marked by dots. Both models exhibit significant dif- 
ferences at nearly all position angles. 

Upper and lower limits for 0,^ are given by the following 
arguments: Since rcR must be located outside the bar, Up must 
not exceed a certain value, depending on the mass distribution 
in the model. This leads to an upper limit for ilp. Estimating the 
lower limit is more complicated since the height of the humps 
in the rotation curves is not exclusively determined by ilp, but 
is also influenced by several other model parameters. 

Examples for this mutual influence of several parameters 
are the cases of models 1 and 2: Numerical N-body-simulations 
of barred galaxies strongly suggest that the CR of intermediate 
type galaxies (e.g. NGC 2336) normally is located at tcr ~ 
1. labor- In models 1 and 2 (Section 9.1), CR is placed at the 
end of the bar at rcR — l.Oobar^ which in fact is the lower 
boundary for tcr- Even for the more realistic model 1 with the 
lower bar mass, the Hll-velocities are still too small. Therefore, 
if we reduce fip of model 1 to a value shifting Corotation to 
rcR = l-la^ar, the velocity humps of the PA = 94°-rotation 
curve increase towards values that are in much better agreement 
with the observations. 

10. The Optimal Model for NGC 2336 

In the last Section we used models with different bar masses 
(1 and 2) and different disk masses (3 and 4) to illustrate the 
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Fig. 12. X2-rotation curves of model 3 (Cdisk ~ 5.0, left col- 
umn) and model 4 (Cdisk — 2.5, right column). Observational 
values are marked by dots. While model 3 produces X2-orbits 
for a neglegibly small central area of NGC 2336 only, much 
more orbits are obtained in model 4. Similarly as in models 1 
and 2, X2-orbits do not reproduce the observed velocities, ei- 
ther, especially for a slit orientation along the disk minor axis 
(PA = 94°). 



qualitative and quantitative changes in the division of the phase 
space and the Hll-kinematics which are caused by the variation 
of Mbar and Mdisk- The effect of Op-variations was described 
as well. But as was discussed in the beginning of Section 9, 
there are many more parameters which may induce changes in 
the predicted kinematics, e.g. the scale lengths of disk, bulge 
and bar, r^, rt,, and a. The results of the studies of changes in 
the latter parameters, as well as the results of the models 1-4, 
can be summarized in the following constraints applying to the 
construction of the final NGC 2336-model: 

1 . Since changing the vertical flattening of the bulge, q, results 
in minor changes in the Hll-rotation curves only, q = 0.4 is 
chosen for the optimal model, which is well in agreement 
with the photometric observations. 

2. The examination of models with different offset angles ijj 
between the bar and the LON shows that a deviation of 
a few degrees (±2°) leads to small changes in the model 
gas kinematics that tend to change the rotation curves in 
the same way as do variations of the bar axis ratio b/a. 
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Since the morphological decomposition of NGC 2336 sug- 
gests ip = 104° as the most probable number and since a 
good fit to the observed velocity field can be achieved with 
this value, ip — 104° will be used for the final model. 
Models with fip being varied support those values that 
place CR at rcn. — l-lobar- Placing rcR outside the bar 
in an interval rcR = 1.0 . . . lAubar is supported by results 



from numerical studies (cf. Athanassoula 1992 and refer- 
ences herein). 

4. A bar-6/a of 0.3 is in good agreement with the kinematic 
requirements and the observations as well. 

5. The bar mass-to-light ratio should not exceed Char — 2.0. 

6. Due to the difficulties in determining the bulge scale length, 
rfc, this parameter cannot be determined exactly from ob- 
servations. The results of the wide range of models exam- 
ined suggests that r^ should lie within rf, = 0.2 . . . 0.3LU 
(2.1...3.1kpc). 

7. The photometrically derived bulge mass {Chuige = 1-0) is 
a reasonably good choice. 

8. The observed disk scale length of r^ = 0.46LC/ (4.76kpc, 
unchanged in all models) agrees well with the requirements 
of the rotation curves, assumed that all other parameters are 
precisely adjusted as well. 

9. With rd unchanged, a small disk mass (i.e. Cdisk < 3.0) 
seems to be a plausible choice when one tries to fit the 
model rotation curves to the observed kinematics in the 
outer disk regions at r > 50" via adjusting the cloud speeds 
on the a; 1 -orbits beyond CR. 

10. In order to adjust the height of the overshootings in the ro- 
tation curves for a slit orientation PA = 94° and to place 
CR at rcR = l.labar, the pattern speed Vtp has to take low 
values < 20km/sec/kpc. 

The final model was optimized by studying 32 parameter vari- 
ations with total number of ~ 9500 closed orbits. Finally, the 
best results were obtained with the parameter set listed in Table 

I 

The total mass of the final model was calibrated according 

to equation 28 and equation 29 with the outer circular velocity 

value from the kinematical observations. We obtain a total mass 

of Mtot = 1-20 • IO^Mq. If we identify for a rough estimate 

the relative mass-to-light-ratios Cdisk and Cbar with the normal 

astrophysical values, we obtain from our optimal model: 



'Mp 



Mbulge = 1.20 • 10 

Mdisk = 9-63 • lu iviQ 
Mhar = 1.13-10i°M, 



10^°M(T 



0- 



(32) 
(33) 
(34) 



Thus we have a bulge/disk ratio of B/D = 0.125 and a 
bar/disk ratio of Bar/D = 0.117. One should keep in mind 
that especially the disk mass contains a certain contribution of 
dark matter which can not be separated within the framework 
of this paper 

Looking at the artificial rotation curves in Fig. Il3| one rec- 
ognizes the excellent agreement between the artificial curves 
and the observed ones. Not only does the optimal model pro- 
duce the correct maximum rotational velocities of the slit ori- 



bar 


tthar [LU/kpc] 
bbar [LU/kpc] 

b/a 

Mbar [Mq] 
^bar 


0.7/7.17 
0.21/2.15 

0.3 

1.13- 10^° 

1.5 


bulge 


7 

rb [LU/kpc] 

<? 

Mbulge [Mq] 
^bulge 


0.9 
0.23/2.36 

0.4 
1.20 ■ 10^" 

1.0 


disk 


ra [LU/kpc] 

Mdisk [Mq] 
Cdisk 


0.46/4.76 

9.63 • 10^" 

2.5 


Qp [km/sec/kpc] 


16.40 



Table 8. Parameter values for the optimal NGC 2336-model. 
All scale lengths are given in model (LU) and physical units 
(kpc), using a distance of 22.9 Mpc for NGC 2336. 



entations PA = 5° and 28°, but also traces the overshootings 
of the PA — 94° -rotation curve correctly. The positions of the 
velocity maxima and minima in the model rotation curves and 
the observed ones are identical. The theoretical rotation curves 
in Fig. |l3| display gaps in the vicinity of rcR, since the orbit in- 
tegrator can not compute orbits near the Lindblad resonances. 
Small discrepancies that remain between observation and 
theory are the central parts of the rotation curves for PA — 5o 
and PA = 28°, where the theoretical slopes rise slightly faster 
than the observed ones. 

11. Discussion 

11.1. Morphological Decomposition 

The morphological decomposition of the J-band image of 
NGC 2336 seems well established, since nearly all parameters 
that were initially derived from the multi-component 2D model 
of the J-band surface brightness distribution appear in the final 
model with the same numerical values. There are justified ex- 
ceptions, e.g. the bulge scale length r^, where the limited spa- 
tial resolution of the images and the small area covered by the 
bulge component prevent us from a more precise determina- 
tion. The moderate increase of the final value of r^ - compared 
to the morphological decomposition - nevertheless agrees well 
with observations. It is also remarkable that the bar-&/a deter- 
mined from observations is nearly the same as that used for the 
optimal model. This is not self-evident, as the inclined position 
of the bar in the sky makes the determination of the deprojected 
bar-6/a difficult. 

Uncertainties of the model parameters become large wher- 
ever the spatial resolution or the area of the A^/i?-detector is 
concerned: The determination of the bulge scale length rt is 



K. Wilke et al.: Mass distribution and kinematics of ttie barred galaxy NGC 2336 



19 



a 



NGC 2336 PA=118' 
\ \ ^ 




NGC 2336 PA=g4^ 
\ \ H 



NGC 3336 PA=3e" 

H h 



— I — 



NGC 2336 PA=5'= 

I I I L 



300 
300 
100 

-100 
-200 
-300 
200 
100 

-100 
-300 
-300 
BOO 
100 

s. -100 
-200 
-300 
300 
^ 100 

i, -100 

-200 

-300 

-100 -50 50 100 

r [arcsec] 

Fig. 13. Observed (points) and artificial rotation curves (filled 
circles) for NGC 2336. Longslit orientations are indicated ac- 
cording to the position angles (PA) of the spectrograph. The 
points at large radii after the gap are produced by the nearly 
circular Xi-orbits outside the CR radius. 



certainly precise to a 50% level only, due to the small bulge area 
(i.e. , the small number of pixels) and the residuals remaining 
after the subtraction of the bulge and disk components. For the 
disk, similar uncertainties are introduced by the limited size of 
the detector chip which omits the outer parts of the disk, lead- 
ing to an uncertainty in rd of w 15%. The luminous masses of 
disk, bulge, and bar can be determined more precisely, espe- 
cially the contribution of the bar is uncertain by 5 — 10% only, 
since the very regular disk structure ensures a reliable decom- 
position. For the photometrically derived values of Mdisk and 
Mbuige, a precision of 15% seems a realistic estimate, due to 
the reasons discussed above. 

The deprojection angles (f>, 6, and tp are reliably deter- 
mined with an error of ±2° at most, since otherwise a mis- 
match between the observed kinematics and that of the optimal 
NGC 2336-model would occur in at least one slit orientation. 
Obviously this is not the case. 

11.2. Analysis of Phase Space 

From the S'OS'-analysis we conclude that the best NGC 2336- 
models show a phase space partition characteristic for barred 
galaxies with a bar contributing roughly 10% to the total po- 



tential: For all realistic parameter sets, an ILR is produced to- 
gether with a:;2 -orbits that dominate the motion in the innermost 
parts of the NGC 2336-models at low i?j-values. 

For energy values Ej > Ei^n, the a;2-orbits disappear, 
only the prograde xi- and the retrograde a:;4-family remain. In 
the vicinity of Ecr, the xi -family dissolves and is replaced 
by higher resonances, e.g. the (4:1)- and the (6:l)-family. The 
(4:l)-resonance ist not very strong, so we conclude that in time 
dependent models the stellar bar of NGC 2336 would be mod- 
erately box-shaped only. The small amount of semi-ergodicity 
leads to the conclusion that a time-dependent model with our 
optimal parameter set would be secularly stable on time scales 
of the bar rotation period. 

Summarizing the results of phase space analysis, the op- 
timal model for NGC 2336 exhibits the typical behaviour of 
disk-dominated barred galaxies with the bar-6/a and ^bar be- 
ing moderately small. 

11.3. Rotation curves 

In our studies, special emphasis was put on the reproduction 
of the gas kinematics in those regions influenced by the bar, 
i.e. the large streaming motions of Hll-clouds orbiting on xi- 
orbits around the center of NGC 2336. However, the overall 
kinematics with the observed rotational velocities in the outer 
regions of the disk had to be described correctly as well. Our 
final model fulfills both requirements, as is obvious from the 
excellent agreement between observed and artificial rotation 
curves in every part of NGC 2336 and for every slit orientation. 

Moreover, the success of adjusting the outer parts of the 
model rotation curves is an evidence for the good deprojection 
performed here. Otherwise the final velocities in the outer disk 
regions could not be fitted simultaneously for all slit directions. 

In general, all characteristic features of all velocity curves 
are traced well by the final model, especially the overshootings 
that can be observed in slit orientations of PA ~ 5°, 28° and 
94°, are produced by the model to full extent. Small differences 
between observed and artificial gas velocities remain in the in- 
nermost region of NGC 2336 and are probably caused by the 
uncertainties in the bulge scale length ri, due to possible errors 
in the morphological decomposition. The precise determina- 
tion of the latter is difficult due to the limited spatial resolution 
of the A^/i?-detector and the small area covered by the bulge 
component. 

From the observational point of view there is no necessity 
for including X2 -orbits in the final sets of closed orbits, because 
a perfect fit to the observed kinematics is already achieved 
by xi -orbits alone. Moreover, the possibility of X2-orbits con- 
tributing significantly to the Hll-kinematics in the inner part 
of the bar are definitely ruled out, since no model with X2- 
orbits can account for the streaming velocities measured at a 
slit orientation of PA = 94°. Though our final model pro- 
duces a moderately strong ILR and therefore numerous X2- 
orbits, this is no contradiction to their uselessness in the case 
of NGC 2336, because our model does not make any statements 
with regard to the occupation number of single orbit families. 
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To correct for the presence of dark matter in the outer re- 
gions of NGC 2336, moderate corrections of the conversion 
factors Cdisk and Cbar had to be applied in order (i) to produce 
the large streaming velocities along ai,ar and (ii) to normal- 
ize the circular xi -velocities in the outer disk regions: Cdisk 
changed from the observationally determined value of 1.0 to 
2.5, while the visible bar mass (Cbar, nir = 1-0) had to be in- 
creased by only 50%. The value Cbar — 1.5 is a further hint at a 
reliable morphological decomposition: For the bar component 
which is not strongly affected by a possible dark matter halo, 
a moderate increase of the NIR-Cbar seems plausible, espe- 
cially when considering possible stellar population changes be- 
tween bulge and bar that are not known in detail. 

Summarizing the results of our kinematical analysis, the 
final NGC 2336 model rotation curves are a very good repre- 
sentation of the observed velocity field. Obviously, there is no 
necessity of including additional hydrodynamical or magneti- 
cal effects in order to explain the observed kinematics. It is also 
remarkable that these results can be achieved based on models 
that neglect spiral arms and an explicit dark halo component 
with independent scale length. 

12. Summary 

In this paper, 2D models of the observed J-band luminosity 
of NGC 2336 were constructed by fitting a disk, a bulge, and 
a bar to the observed surface brightness distribution. The re- 
sulting model was deprojected and converted to an underlying 
mass distribution. The total 2D potential, ^ total, was obtained 
by numerical expansions of the potential of the single compo- 
nents and coadding them. For the examination of test particle 
motions in the potential of NGC 2336, we used a numerical 
orbit integrator, in which ^totai and its derivatives are imple- 
mented. The time-dependent motions of single Hll-clouds are 
traced by a grid-based integration scheme, in which particles 
move on arbitrary orbits in the stationary potential of the rotat- 
ing bar. 

The resulting partition of the phase space was examined by 
analysing the appropriate cuts through phase space at a given 
energy, the Poincare Surfaces of Section. They provided us 
with the start values - coordinates and velocities - of the main 
orbit families constituting the bar. Given those initial values, 
complete sets of closed orbits in a certain energy range were 
computed and projected to the sky. Using a virtual longslit in 
various orientations, artificial rotation curves were constructed 
and compared to the observed kinematics of the Hll-gas. 

In a further step, the effects of parameter variation on the 
division of phase space and the predicted Hll-rotation curves 
of numerous NGC 2336-models were examined. The effects of 
varied Mbar- and M^igfe-values were shown here in greater de- 
tail, while the effects of all other free parameters being changed 
and studied are described qualitatively. The models led to nar- 
row limits for the parameter values of an optimal model which 
should (i) reproduce the observed Hll-kinematics directly and 
(ii) use the parameter values directly obtained from the mor- 
phological decomposition procedure wherever possible. 



From the final NGC2336-model, the following results are 
obtained: 

1. The model parameter values obtained from observations 
and those needed for an optimal fit of the Hll-kinematics 
are identical (with justified exceptions), therefore the aim 
of constructing a consistent NGC 2336-model is achieved. 

2. The overall quality of the morphological decomposition 
model is sufficiently good to produce a realistic mass model 
of NGC 2336. Nevertheless, to check whether the kine- 
matically supported scale length values of bulge and disk 
are fully consistent with the more uncertain morphologi- 
cal ones, we would need Af/_R-data with (i) a better spatial 
resolution and (ii) larger areas covered by the detector chip. 
With the small pixel number and the comparably large pixel 
size of the MAG/C-detector, the determination of r^ and 
especially rf, remains crucial. 

3. The phase space of all acceptable models, especially the 
optimal one, exhibits a structure which is typical for a disk- 
dominated galaxy with a bar roughly contributing 10% to 
the total potential, i.e. small amounts of semi-ergodicity as 
well as the presence of the xi-, X2- and a;4 -family of peri- 
odic orbits. 

4. a:;2-orbits can definitely be excluded when one tries to fit 
the observed velocities within the bar, since the streaming 
velocities produced by them are by far too low. 

5. The optimal model favours a low pattern speed {Vtp w 
16km/sec/kpc), a small bar axis ratio {h/a — 0.3) and - 
compared with the luminous matter - moderately increased 
bar and disk masses due to the presence of dark matter 

{Mbar = 1.13 • lOlOMo, Md^sk = 9.63 • WMq). 

6. An excellent agreement between the observed and the ar- 
tificial rotation curves is achieved, the optimal NGC 2336- 
model describes correctly the large streaming motions of 
the Hll-clouds along the bar as well as the rotation veloci- 
ties in the outer regions of the disk. 

7. In general, the case of NGC 2336 shows that the peculiar 
Hll-kinematics of a strongly barred galaxy can be com- 
pletely explained by consistent stationary models which are 
constructed from observationally derived parameter sets. 

Acknowledgements. The authors want to thank R. Bender (Univer- 
sitatssternwarte Miinchen) for providing the software for the Fourier 
quotient correlation analysis of the stellar spectra, FCQ4. We also ac- 
knowledge the help of J. Heidt with the NIR data reduction. This 
work has been supported by the DFG Sonderforschungsbereich 328. 

References 

Athanassoula, E., 1992, MNRAS 259, 345 

Bender, R., 1990, A&A 229, 441 

Binney, J., Tremaine, S., 1987, Galactic Dynamics, Princeton Univer- 
sity Press 

Combes, F., 1994, in: The Formation of Galaxies, Munoz-Tunon C. 
(ed.),5th Canary Island Winter School 

Contopoulos, G., Papayannopoulos, Th., 1980, A&A 92, 33 

Dehnen, W., 1993, MNRAS 265, 250 

de Vaucouleurs, G., 1948, ARA&A 11, 247 



K. Wilke et al.: Mass distribution and kinematics of ttie barred galaxy NGC 2336 21 



Ferrers, N.M., 1877, QJPAM 14, 1 

Frogel, J.A., Quillen, A.C., Pogge, R.W., 1996, in: New Extragalactic 
Perspectives in New South Africa, D.Block (ed.), Kluwer, Dor- 
drecht, 65 

Goldstein, H., 1981, Classical Mechanics, Addison- Wesley Publish- 
ing Company, 145 

Hernquist, L., 1990, ApJ 356, 359 

Jaffe, W., 1983, MNRAS 202, 995 

Martin, P, 1995, AJ 109, 2428 

Matthias, M., 1993, diploma thesis. University of Heidelberg 

Salucci, P, Ashman, K., Persic, M., 1991, ApJ 379, 89 

Sellwood, J.A., Wilkinson, A., 1992, Reports on Progress in Physics 

Toomre, A., 1962, ApJ 138, 385 

Tully, R., B., 1988, Nearby Galaxies Catalogue, Cambridge University 
Press 

van Moorsel, G.A., 1983, A&AS 54, 19 



